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1 Introduction 

This paper presents a search for physics beyond the standard model (SM) in events containing a 
pair of opposite-sign same-flavor (SF) electrons or muons, jets, and an imbalance in transverse 
momentum. The analysis is based on a sample of proton-proton (pp) collisions collected at 
a center-of-mass energy of 8TeV with the CMS detector 0] at the CERN LHC in 2012 and 
corresponds to an integrated luminosity of 19.4 fb _1 . 

The invariant mass distribution of the two-lepton system can exhibit an excess that increases 
with the dilepton mass, followed by a sharp decrease and thus an "edge", if the two leptons 
originate from the decay of an on-shell heavy neutral particle. This kind of signature is fairly 
generic for models of physics beyond the SM (BSM), assumes an isotropic decay, and is purely 
kinematic in origin. In models of supersymmetry (SUSY) f2j, an edge with a triangular shape 
is expected in the cascade process X 2 —> li —> y < l ) T f U 0, where X 2 ar >d x\ are respectively the 
next-to-lightest and lightest neutralino, with £ a selectron or smuon, the SUSY partners of an 
electron or muon. Alternatively, the X 2 can undergo a three-body decay to X^£ + £~ through a 
virtual Z* boson, also yielding an edge in the dilepton mass spectrum but with a more rounded 
shape. Another possibility is the decay of a ^ to an on-shell Z boson, Xi ~t z'j’Z. This latter 
process does not produce an edge but rather a dilepton mass peak near 91 GeV. These processes 
arise as a consequence of the gauge-coupling structure of SUSY and are a characteristic feature 
of SUSY decay chains. Their relative importance depends on the SUSY mass hierarchy and is 
thus model dependent. 

This search is therefore motivated by the possible existence of the fairly generic signal shape of 
an edge, or of a peak at the Z boson mass, that would be visible in the invariant mass distri¬ 
bution of the two leptons. The position of the edge would give an indication of the unknown 
BSM mass hierarchy. The dilepton invariant mass provides a search variable that is unaffected 
by uncertainties in the jet energy scale and resolution, and the search for an edge is therefore 
complementary to searches based solely on hadronic quantities. 

The CMS Collaboration previously presented two searches for BSM physics based on events 
with an opposite-sign SF lepton pair: a search for an edge in the dilepton mass spectrum out¬ 
side the Z boson mass region 0), and a search for events containing on-shell Z bosons |5|. Both 
these studies were conducted using the 7 TeV CMS data sample collected in 2011. The present 
study updates and combines these two analyses, using the 8 TeV data sample. Searches for 
SUSY in opposite-sign dilepton final states, but targeting different production and/or decay 
mechanisms, are presented by the ATLAS Collaboration in Refs. ©El- 

A brief description of the CMS detector is given in Section [2] Signal models studied in this 
analysis are described in Section [3] Section |4] outlines the event selection and simulation. The 
background estimation methods are presented in Section |5j A fitting procedure used to search 
for an edge is described in Section [6] The results of the search are presented in Section [7] 
Systematic uncertainties associated with the signal modeling are discussed in Section[8]and the 
results of the search are interpreted in the context of the signal models in Section[9] A summary 
is presented in Section [To| 


2 Detector and trigger 

The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter 
that produces an axial magnetic field of 3.8 T. A silicon pixel and strip tracker, a lead tungstate 
crystal electromagnetic calorimeter, and a brass/plastic-scintillator hadron calorimeter are po- 
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sitioned within the field volume. Iron and quartz-fiber hadron calorimeters are located outside 
the magnetic field volume, within each endcap region of the detector. Muons are measured 
using gas-ionization detectors embedded in the steel flux-return yoke outside of the solenoid. 
A detailed description of the CMS detector, its coordinate system, and the main kinematic vari¬ 
ables used in the analysis can be found in Ref. [Tj. 

Events must satisfy at least one of a set of ee, pp, and eu triggers. The ee and pp triggers collect 
signal candidate data while the ep trigger collects data used in the background-determination 
procedure, as described below. These triggers require at least one electron or muon with trans¬ 
verse momentum pj > 17GeV, and another with pj > 8GeV. Their efficiencies after event 
selection (>90%) are measured in data and are accounted for in the analysis. The efficiencies of 
the ee, pp, and ep triggers differ by only a few percent in the kinematic range of this search. 


3 Signal scenarios 

Two classes of signal events are considered, as explained below. Both classes are implemented 
in the framework of simplified models 03, in which only the targeted production and decay 
schemes are examined, with all non-participating BSM particles assumed to be too heavy to be 
relevant. 

The first class of signal events targets the production of an edge in the invariant mass spectrum 
of opposite-sign SF lepton pairs, as expected from the correlated production of these leptons 
in cascade decays. This class of scenarios is based on the production of a bottom squark- 
antisquark pair. Each bottom squark b decays to a bottom quark b and the Xi neutralino. 
Two specific possibilities are considered. In the first scenario (Fig. [l] left), the decays to an 
off-shell Z boson Z* and the y', 1 neutralino, where the ^ is a stable, weakly interacting, lightest 
SUSY particle (LSP). The Z* boson decays according to its SM branching fractions, sometimes 
producing a charged lepton pair £ + £~ (£ = e, p). The mass difference between the x { { and x'\ / 
which determines the location of the edge, is fixed to 70 GeV. This scenario is referred to as the 
"fixed-edge" scenario. In the second scenario (Fig. [l] right), the yfj decays to an on- or off-shell 
Z boson and the x'l LSP or according to x% —> ££, with a 50% probability for each decay. The 
slepton £, i.e., the SUSY partner of the lepton, then decays according to l —> £^\. The consid¬ 
ered sleptons are mass-degenerate selectrons and smuons. The mass of the slepton is chosen to 
lie halfway between the masses of the two neutralinos: nij = m ~o + 0.5(m-o — m-o). The mass 
of the x_ ( \ is set to 100 GeV, with the position of the edge a free parameter in a scan of the mass 
spectrum. This scenario is referred to as the "slepton-edge" scenario. 

The second class of signal events targets the production of an opposite-sign SF lepton pair from 
the decay of an on-shell Z boson. This class of scenarios, illustrated in Fig. |2j is based on gluino 
pair production in the context of gauge mediated supersymmetry breaking (GMSB) models @- 
□n. Each gluino decays to a quark-antiquark pair and the neutralino. The x^ decays to 
an on-shell Z boson and a stable, massless, weakly interacting gravitino LSP. We refer to this 
scenario as the "GMSB" scenario. 

The production of squark and gluino pairs is simulated with the MadGraph 5.1.3.30 IH21H3 I 
Monte Carlo (MC) leading-order event generator, including up to two additional partons at 
the matrix element level. The decays of the squarks, gluinos, and other particles are simu¬ 
lated with the PYTHIA 6.4.22 ffT4j| event generator. The MadGraph events are subsequently 
processed with the PYTHIA program to generate parton showers and account for hadroniza- 
tion. The decay of the Z boson is handled in PYTHIA. In all the aforementioned scenarios 
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Figure 1: Event diagrams for the (left) "fixed-edge", and (right) "slepton-edge" scenarios, with 
b a bottom squark, X 2 the second lightest neutralino, x { \ a massive neutralino LSP, and i an 
electron- or muon-type slepton. For the slepton-edge scenario, the Z boson can be either on- or 
off-shell, while for the fixed-edge scenario it is off-shell. 
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Figure 2: Event diagram for the "GMSB" scenario, with g a gluino, X\ the lightest neutralino, 
and G a massless gravitino LSP. 

(fixed-edge, slepton-edge, GMSB), the Z boson decays according to its SM branching frac¬ 
tions. To reduce computational requirements, the detector response is simulated using the 
CMS fast simulation fl5| . Differences in the lepton reconstruction and identification efficien¬ 
cies between the fast and a "full" simulation, where the full simulation is based on process¬ 
ing through the GEANT4 Itl6l programs, are corrected using scale factors. The expected signal 
event rates are normalized to cross sections calculated at the next-to-leading order (NLO) in the 
strong coupling constant, including the resummation of soft gluon emission at next-to-leading- 
logarithmic (NLO+NLL) accuracy lfl7H22ll . 


4 Event selection, reconstruction, and search strategy 

We select events with an oppositely charged lepton pair (e 1 e e ± p T , or p + p~). The leptons are 
required to have pj > 20GeV and >/| < 2.4, where rj is the pseudorapidity [T). In events with 
more than two selected leptons, we choose the two oppositely charged leptons with highest pj. 
The background estimation techniques employed in this analysis rely, in part, on a symmetry 
between lepton pairs with the same flavor compared to those with opposite flavor (OF), where 
OF refers to e~p combinations. It is therefore desirable that the reconstruction efficiencies of 
electrons l23l and muons |[24l be as similar as possible. For this reason, we exclude leptons 
in the intervals 1.4 < \rj\ < 1.6 between the barrel and endcap regions of the detector (U, 
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where the electron and muon reconstruction efficiencies differ significantly. For the SF signal 
candidate sample, only events with an e + e _ or p' p pair are used. 

Leptons produced in the decays of low-mass particles, such as hadrons containing b and c 
quarks, almost always lie in or near jets. The background from these low-mass processes can 
be suppressed by requiring the leptons to be isolated in space from other particles. A cone of 
radius A R = \/(A p) 2 + (A cp) 2 = 0.3 is constructed around the lepton momentum direction, 
where (p is the azimuthal angle. The lepton relative isolation is then quantified by the scalar 
px sum of all particle candidates within this cone, excluding the lepton, divided by the lepton 
pj. The resulting quantity is required to be less than 0.15. The sum includes a correction 
to the total energy to account for contributions from extraneous pp interactions within the 
same or a nearby bunch crossing (pileup). For electrons, the pileup correction is based on the 
effective area method Il25l . while for muons it is based on the number of charged hadrons not 
associated with the primary vertex. The performance of the electron and muon isolation criteria 
is discussed in Refs. Il23l[24 |. 

The primary vertex is taken to be the reconstructed vertex with the largest p\ sum of associated 
tracks. Leptons with impact parameters larger than 0.2 mm in the transverse plane or 1 mm 
along the beam direction are rejected. As an additional requirement, the two selected leptons 
must be separated by A R > 0.3 to avoid systematic effects that arise for isolation requirements 
in very collinear topologies. 

A particle-flow (PF) technique Il26l is used to reconstruct jets and missing transverse momen¬ 
tum. All objects reconstructed with the PF algorithm serve as input for jet reconstruction, based 
on the anti-kx clustering algorithm H27l with a distance parameter of 0.5, as implemented in the 
FASTjET package [[25] [28J. We apply pj- and //-dependent corrections to account for residual 
effects of nonuniform detector response. The contribution to the jet energy from pileup is es¬ 
timated on an event-by-event basis using the jet area method described in Ref. f29i . and is 
subtracted from the overall jet pj. Jets are required to have pj > 40GeV, |//| < 3.0, and to be 
separated by A R > 0.4 from all selected leptons. The missing transverse momentum pj" lss is 
defined as the projection on the plane perpendicular to the beam axis of the negative vector 
sum of the momenta of all reconstructed PF objects in an event. The magnitude of p™ iss is 
referred to as E™ ISS . The E™ 1SS distributions of events in the SF and OF samples for dilepton 
invariant mass m^> 20GeV and number of jets Nj e t s > 2 are shown in Fig. |3j 

The event selection criteria are motivated by the expectation that BSM signal events, involving 
the production of new heavy particles, generally have larger jet multiplicity and Ej' ss than 
background events, which primarily arise from top quark-antiquark (tt) and Drell-Yan (DY) 
processes. The large value of Ej' ss expected in signal events is due to the weakly interacting 
LSP particles, which escape without detection. 

In the search for an edge, we define two signal regions: either Nj ets > 2 and E™ ISS > 150 GeV, or 
Njets A 3 and E™ ISS > 100 GeV. For both regions, we separately consider events in which both 
leptons satisfy |//i ep | <1.4 ("central" signal region) and events in which at least one lepton 
satisfies 1.6 < |//i e p| < 2.4 ("forward" signal region). The motivation for the distinction be¬ 
tween the central and forward regions is that for BSM production through the decay of heavy 
resonances, the final-state leptons and jets are expected to be more centrally distributed than 
for the SM backgrounds. Two methods are used to search for an edge signature. In the first 
method, a search for an edge is performed in the range 20 < /////. < 300 GeV by fitting the signal 
and background hypotheses to data, as described in Section[6j In the second method, based on 
a direct comparison of event counts, with no assumption about the shapes of the signal and 
background distributions, we select three regions, 20 < niff < 70 GeV, 81 < < 101 GeV, 
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CMS 19.4 fb" 1 (8 TeV) CMS 19.4 fb" 1 (8 TeV) 



Ej iss [GeV] E™ 33 [GeV] 


Figure 3: The Ej' ss distributions of events in the SF (left) and OF (right) samples for nifg > 
20 GeV and Nj ets > 2 in comparison with predictions for the SM background from the MC 
generators described in Section[4] In the ratio panel below each plot, the error bars on the black 
points show the statistical uncertainties of the data and MC samples, while the shaded band 
indicates the MC statistical and systematic uncertainties added in quadrature. The rightmost 
bins contain the overflow. 


and ni(( > 120 GeV, denoted the "low-mass", "on-Z", and "high-mass" regions, respectively. 
For this "counting experiment", the integrated yield in each region is compared to the corre¬ 
sponding background prediction. 

In the search for BSM events with an on-shell Z boson, we perform a dedicated counting exper¬ 
iment in the region 81 < nifj < 101 GeV, based on the distribution of Ey' ss . For this study, two 
inclusive bins in the number of jets are defined: Nj e t s > 2 and Nj e t s > 3. Events are examined 
in exclusive bins of E™ 1SS as described in Section [zj 

While the main SM backgrounds are estimated using data control samples, simulated MC 
events are used to evaluate smaller sources of background. The simulation is also used to esti¬ 
mate uncertainties. Simulated samples of DY + jets, tt, W + jets, VV, and ttV (V = W, Z) events 
are generated with the MadGraph 5.1.3.30 event generator. The DY + jets sample considers 
events with dilepton invariant masses as low as 10 GeV, as well as decays to the tt channel. 
The matrix element calculation performed with MadGraph 5.1.3.30 is then interfaced to the 
PYTHIA 6.4.22 program for the description of parton showering and hadronization, analogous 
to the procedure used for the signal samples. The detector response in these samples is simu¬ 
lated with a Geant4 model fl6l of the CMS detector. The simulated events are reconstructed 
and analyzed with the same software used to process data. In the simulation, multiple pp in¬ 
teractions are superimposed on the hard collision, and the simulated samples are reweighted 
to reflect the beam conditions, taking into account the total inelastic pp cross section at the 
LHC. Corrections are applied to account for the differences between simulation and data in 
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5 Background estimates 


the trigger and reconstruction efficiencies. The simulated sample yields are normalized to an 
integrated luminosity of 19.4 fb 1 using NLO inclusive cross sections, except for the DY + jets 
and W + jets samples, where next-to-next-to-leading order calculations 13011 are used. 


5 Background estimates 

The principal SM backgrounds are divided into two categories. Backgrounds that produce 
OF pairs (e + p~,e~p + ) as often as SF pairs (e 1 e - ,p 1 p ) are referred to as flavor-symmetric 
(FS) backgrounds. This category is dominated by tt processes. Drell-Yan events form the 
second principal background category. The FS background estimate accounts also for WW, 
Z/ 7 *(—> tt), and tW single-top quark production, as well as for backgrounds due to leptons 
from hadron decays and from hadrons misidentified as leptons. Contributions from tt + X, 
with X a W, Z, or Higgs boson, have been studied in the simulation and found to be negligible. 

The missing transverse momentum in DY + jets events arises primarily from jet energy reso¬ 
lution and reconstruction effects. Contributions from SM WZ and ZZ processes, which might 
include genuine E™ 1SS , are incorporated into the DY + jets background estimates. 

5.1 Flavor-symmetric backgrounds 

The contribution of FS background events to the signal regions is determined using a control 
sample defined by OF events that satisfy the full event selection. The OF yields are multiplied 
by a factor Rg F /or to account for efficiency differences in the selection of dilepton pairs in the 
SF and OF samples. Two methods are used to evaluate the Rg F /of factor, as explained below. 
We also determine corresponding factors R ee /OF and R ;(? ,/of for the individual SF terms. 

The first method, referred to as the "factorization" method, uses a data control sample to de¬ 
termine the SF-to-OF ratio of reconstruction efficiencies, and a second data control sample to 
determine the corresponding ratio of trigger efficiencies. The reconstruction efficiency ratio is 
evaluated using a large DY + jets data control sample selected by requiring 60 < mgg < 120 GeV, 
A/jets > 2, and Ej' ss < 50 GeV. The measured quantity is r }ie = ^/N^^/N ee , which represents a 
ratio of the muon-to-electron reconstruction efficiency uncorrected for trigger efficiencies. The 
extrapolation of r f/e into the central (forward) signal region is studied using data. A 10% (20%) 
systematic uncertainty is assigned to account for the observed dependencies. The ratio Rj of 
the trigger efficiencies in SF and OF events is determined using a data control sample selected 
with a trigger based on Hj, which is the scalar sum of jet fj values for jets with pj > 40 GeV. 
The sample is selected by requiring Hj > 200 GeV, mgg > 20 GeV, and excluding events with 
Njets = 2 and E™ ISS > 100 GeV. The latter two requirements ensure that this control sample is 
uncorrelated with the control sample used to determine the reconstruction efficiencies. A 5% 
uncertainty is assigned to the efficiency of each dilepton trigger (e + e _ , p + p~, e ± p : ) to account 
for the observed dependence of these efficiencies on mgg, E™ ISS , and Nj e t s - The full correction to 

the OF event rate is R SF /of = 0.5(r,, e + r~^)R T , with R T = \Je^e^/e^. 

The second method, referred to as the "control-region" method, directly measures the Rg F /of 
factors from the ratio of the e + e~, p + p~, or combined SF yields with the OF yields in a tt- 
dominated control region with Nj e t s > 2, 100 < E™ 1SS < 150 GeV, and 20 < mgg < 70 GeV. 
Differences between the factors obtained in this control region and in the signal regions are 
studied with the tt simulation. No difference is observed within the MC statistical uncertainty, 
which is used to define the corresponding systematic uncertainty. The size of this uncertainty 
is found to be 2% (3-4%) in the central (forward) region. 





5.2 SM Drell-Yan background 
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As the two methods rely on uncorrelated control samples, and the results agree within their 
uncertainties, the final results for Rsf/of are obtained by taking the weighted average of the 
results from the two methods, propagating the uncertainties. The results are shown in Table [l] 
The FS background is then given by the number of events in the OF control samples multiplied 
by the corresponding Rg F /q F factor. It is seen that the values of the Rgp /of ratios are consistent 
with unity within the uncertainties. We find that the variation of Rsf/of for increasing m^, 
A/jgts, and E™' ss lies within the assigned uncertainty We thus use the same value of Rgp /of for 
all signal regions. 

Table 1: Results for Rsf/of i n the signal regions. The results of the two methods are shown 
with statistical and systematic uncertainties, while the uncertainties for the combined values 
are a combination of the statistical and systematic terms. The values of r 1(C and Rj listed for 
the control-region method are not used directly in the analysis and are listed for purposes of 
comparison only 



Central 

Forward 

Factorization method 

Rsf/of 

1.03 ±0.01 ±0.06 

1.11 ±0.04 ±0.08 

Ree/OF 

0.47 ± 0.01 ± 0.06 

0.46 ±0.02 ±0.10 

R}i}i /OF 

0.56 ± 0.01 ± 0.07 

0.65 ±0.03 ±0.14 

r }ie 

1.09 ±0.00 ±0.11 

1.18 ±0.00 ±0.24 

Rj 

1.03 ±0.01 ±0.06 

1.10 ±0.04 ±0.07 

Control-region method 

Rsf/of 

0.99 ± 0.05 ± 0.02 

1.11 ±0.11 ±0.03 

Ree/OF 

0.44 ± 0.03 ± 0.01 

0.49 ± 0.06 ± 0.02 

Rf(/(/OF 

0.55 ± 0.03 ± 0.01 

0.62 ±0.07 ±0.02 

T}ie 

1.12 ± 0.04 (stat) 

1.12 ±0.08 (stat) 

Rj 

0.98 ± 0.05 (stat) 

1.11 ± 0.11 (stat) 


Combined 


Rsf/of 

1.00 ±0.04 

1.11 ±0.07 

Ree/OF 

0.45 ± 0.03 

0.48 ± 0.05 

Rf</(/OF 

0.55 ± 0.03 

0.63 ± 0.07 


5.2 SM Drell-Yan background 

In the search for an edge based on a fit, the DY ± jets background is determined as described 
in Section[ 6 j For the counting experiment method, the DY background is determined using the 
jet-Z balance (JZB) and E™ 1SS -template methods QQ, as described in this section. 

The JZB is a measure of the imbalance between the pj of the Z/ 7 * boson and the pj of the 
recoiling hadronic system in DY + jets. The JZB is defined as the scalar difference between the 
pj of the Z/ 7 * and the net pj of the recoiling hadronic system. Standard model DY ± jets events 
equally populate negative and positive values of JZB, because non-zero JZB in these events 
arises from jet energy resolution effects, whereas in BSM and tt events, which contain genuine 
Ej' ss , JZB can be very asymmetric towards positive values because of the correlated production 
of the lepton pair and the undetected particles. Events with negative values of JZB mainly 
arise from DY + jets processes, with a small contribution from tt production. The number of tt 
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events that contribute negative JZB values is determined using the corresponding results for OF 
events. The tt contribution is then subtracted from the number of negative JZB events in the SF 
sample to estimate the DY+jets background. Uncertainties arising from imperfect knowledge 
of the Rgp /of factor (Table ml when subtracting the tt contribution are propagated to the final 
DY + jets estimate. This method also accounts for processes with DY—> + X, where X 

denotes other particles that might be present in the final state. The systematic uncertainty in 
the assumption that DY + jets events equally populate positive and negative values of JZB is 
evaluated in simulation by comparing the E™ 1SS distributions of events with JZB < 0 and >0. 
A systematic uncertainty of 20% is assigned to account for possible differences, dominated by 
the limited statistical precision of the MC sample. 

The £™ ss -template method relies on a data control sample consisting of events with photons 
and jets to evaluate the DY + jets background in a high E™ ISS signal region. For both DY + jets 
and 7 + jets events, large values of E™ ISS arise from the mismeasurement of jet pj values. The 
Ey ISS distribution for DY + jets events can thus be evaluated using 7 + jets events selected with 
similar kinematic requirements. After selection, the Nj e t s , Hj, and boson pj distributions of the 
7 + jets sample are reweighted to match those of the DY + jets sample. The systematic uncer¬ 
tainty is determined in simulation by applying this reweighting to a 7 + jets MC sample and 
comparing the reweighted E™ ISS spectrum to that in a DY + jets MC sample. The uncertainty is 
taken as the larger of the difference between the two samples or the MC statistical uncertainty. 
The assigned uncertainty is 4% for E™ ISS < 60 GeV, 15% for 60 < Ey" ss < 200 GeV, 34% for 
200 < Ey ISS < 300 GeV, and 100% for Ey ISS > 300 GeV. The uncertainty in the last two Ey ISS 
bins is mostly due to the MC statistical uncertainty. Uncertainties are also assigned to account 
for the difference in the number of pileup interactions in the triggered DY + jets and 7 + jets 
events (1-3%, increasing with E™ 1SS ) and the purity of the photon selection (1-5%, increasing 
with E™ iss ). 

With the E™ 1SS -template method, the backgrounds from WZ, ZZ, and other rare SM processes 
are estimated using simulation. An uncertainty of 50% is assigned to the WZ and ZZ back¬ 
grounds based on comparisons with data in orthogonal control samples selected by requiring 
A/jets > 2 and either exactly three leptons (WZ control sample) or exactly four leptons (ZZ con¬ 
trol sample). For the other rare backgrounds, which include events with ttZ, ZZZ, ZZW, and 
ZWW production, an uncertainty of 50% is similarly assigned. 

The two methods are used to estimate the yield of SM DY events in the on-Z region, with all 
other selection criteria the same as for signal events. Since the two methods use uncorrelated 
samples to describe the high-Ey ,ss tail of SM DY + jets production, the DY background estimate 
in the on-Z region is taken to be the weighted average of the two estimates. The individual 
results are consistent with the weighted average within their uncertainties. 

This estimate of the DY + jets background is extrapolated outside on-Z region using the ratio 
Rout/in- h is measured in data as the event yield outside the on-Z region divided by the yield in 
the on-Z region, for the dilepton invariant mass distribution in SF events with Ey 1SS < 50 GeV 
and Njets > 2, after subtraction of the FS backgrounds. This distribution is almost entirely com¬ 
posed of lepton pairs from DY processes, making it suitable for the determination of R ou t/m- A 
systematic uncertainty of 25% is assigned to R out /i n to account for possible biases introduced 
by the different selection criteria used for the signal and control regions, and for differences in 
the value of R out /i n between electrons and muons. 
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6 Kinematic fit 


The search for an edge based on the fit method is performed using a simultaneous extended 
unbinned maximum likelihood fit to the dilepton mass distributions of e ! e , ft' ft ", and e^ft 
events. The likelihood model contains three components: (a) an FS background component, 
(b) a DY background component, and (c) a signal component. 

The FS background is described using a model with three regions: a low-mass region model¬ 
ing the rising distribution shaped by lepton acceptance requirements, a transition region, and 
a high-mass region in which the distribution falls exponentially. The probability density func¬ 
tions describing the three regions are: 


'Pfs^u) = 


, 'P FS ,i(mee) = cim a u 

'PfsAwu) = E?=i c 2,i m'w 1 

.'PFS^( m u) = C3e~P mu 


if 20GeV < met < m u> 
if m$ < m u < m$, 
if mff < m u < 300 GeV, 


( 1 ) 


where ni^ and mffi define the boundaries between the regions. Because of the requirement 
that the function and its derivative both be continuous, the FS model is left with five indepen¬ 
dent parameters plus the normalization. 

The DY background is modeled with the sum of an exponential function, which describes the 
low-mass rise, and a Breit-Wigner function with a mean and width set to the nominal Z boson 
values |3ll , which accounts for the on-Z lineshape. To account for the experimental resolu¬ 
tion, the Breit-Wigner function is convolved with a double-sided Crystal-Ball H32l function 

'Pdscb{™m) : 


where 


'Pdscb^u) = 


Ai(Bx- 


m ll—l l CB \—th 


a CB ‘ 

p Vn i 

exp V 2 0 *B 
a 2( b 2 + mb; 


it 

it 


-n 2 


C'CB 


< -«1, 


- ai <"™< a2 . 


mg-FCB 

CCB 


> 0t 2 , 



and Bi = —r 


Oti 


The full model for the on-Z DY lineshape is thus: 


( 2 ) 


(3) 


"P DY, on-Z (jM-ii) = J 'PoSCB^n^BW^U ~ m') dm'. 


(4) 


The signal component is described by a triangular shape, convolved with a Gaussian distribu¬ 
tion to account for the experimental resolution: 


^s( m a) = 


1 

yjlmia 


edge 

tl 


V exp 


(>»« - yf 

2cr « 


d y- 


(5) 


As a preliminary step, a fit is performed separately for electrons and muons in the DY-enriched 
control region (the same control region as described for R ou t/i n in Section [T2| to determine 
the shape of backgrounds containing a Z boson. In systematic studies of the fit, the DY shape 
parameter that models the exponentially falling spectrum of the virtual photon is varied by 
±25%. The total effect on the fitted signal yield is found to be negligible. The parameters of 
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7 Results 


the DY shape are then fixed and only the normalizations of these backgrounds are free param¬ 
eters in the fit. The nominal fit is applied simultaneously in the signal region to the dilepton 
invariant mass distributions in the e + e _ , }i 1 pi , and samples. Therefore the model for 

the FS background is the same for the SF and OF events. The Rgp /of factors, in the central and 
forward regions, are treated as nuisance parameters, parametrized by Gaussian distributions 
with a mean value and standard deviation given by the values of Rgp /of and their uncertain¬ 
ties (Table [!}. The fit is carried out in the central and forward regions simultaneously, with 
the position of the edge as the only shared parameter. Therefore the signal model has three 
free parameters: the fitted signal yield in the central and forward regions separately, and the 
position of the edge. 


7 Results 

The dilepton mass distributions and the results of the fit in the central and forward signal 
regions are shown in Fig. [4] Table [2] presents a summary of the fit results. A signal yield 
of 126 ± 41 (22 ± 20) events is obtained when evaluating the signal hypothesis in the central 
(forward) region, with an edge located at 78.7 ±1.4 GeV. The p-value, evaluated using —2 In Q, 
is 0.009, where Q denotes the ratio of the fitted likelihood value for the signal-plus-background 
hypothesis to the background-only hypothesis for the case where the edge position is fixed to 
the observed value. This p-value is interpreted as the one-sided tail probability of a Gaussian 
distribution and corresponds to an excess in the observed number of events compared to the 
SM background estimate of 2.4 standard deviations. 

As cross-checks, we tested alternative shapes for the FS background, specifically the sum of 
three Gaussian distributions, and binned and smoothed histograms taken from the OF events. 
In all cases, the results were found to be consistent with the nominal results. 

Table 2: Results of the unbinned maximum likelihood fit for event yields in the signal regions. 
The quoted uncertainties are calculated using the MINOS B33l program and account for both 
statistical and systematic sources. 



Central 

Forward 

Drell-Yan 

158 ± 23 

71 ±15 

OF yield 

2270 ± 44 

745 ± 25 

Rsf/of 

1.03 ± 0.03 

1.02 ±0.05 

Signal events 

126 ±41 

22 ±20 

m edge 

m u 

78.7 ±1.4 GeV 

Local significance 

2.4 j 


Besides the fit described in Section [6j we perform a counting experiment in the mass windows 
20 < mu < 70 GeV, 81 < < 101 GeV, and > 120 GeV, with no assumption about a 

particular signal shape, as mentioned in Section [4] Figure [5] shows the invariant mass distri¬ 
butions for the signal candidate sample and the estimated background. For the background 
prediction, the OF yield in the signal mass window is multiplied by the Rsf/of factor, and the 
background prediction for backgrounds containing a Z boson in the on-Z region by the R 0 ut/in 
factor, as described in Sections [5T] and [±2] 

The results are summarized in Table [3] The significance of the excess in the observed number 
of events with respect to the estimated number of SM background events is evaluated using 
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3 

CL 


m„ [GeV] 



m„ [GeV] 



Figure 4: Fit results for the signal-plus-background hypothesis in comparison with the mea¬ 
sured dilepton mass distributions, in the central (top) and forward (bottom) regions, projected 
on the same-flavor (left) and opposite-flavor (right) event samples. The combined fit shape 
is shown as a blue, solid line. The individual fit components are indicated by dashed lines. 
The flavor-symmetric (FS) background is displayed with a black dashed line. The Drell-Yan 
(DY) background is displayed with a red dashed line. The extracted signal component is 
displayed with a green dashed line. The lower plots show the pull distributions, defined as 

(^data ^fit)/Gdata- 
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7 Results 


a profile likelihood asymptotic approximation It34l . The local significance of the excess in the 
central low-mass region, where the largest deviation is found, is 2.6 standard deviations. Note 
that the signal regions were defined before the data sample was examined, and that the low- 
mass region (20 < < 70 GeV) does not include events between 70 GeV and the best-fit value 

for the location of the edge {m^ = 78.7 GeV). The flavor of the i + i~ pair was studied in the 
counting experiment. Within the statistical uncertainty and accounting for differences in the 
reconstruction efficiencies, the electron and the muon channels are found to contribute evenly 
to the excess. Further studies of the excess in the low-mass region do not yield evidence for a 
neglected systematic term. The excess is observed predominantly in events with at least one 
identified bottom quark jet (b jet) and diminishes if a veto on the presence of a b jet is applied. 
To identify b jets, we use the CSV algorithm at the medium working point l35l . 

Figure[6]presents a comparison of the data and the SM simulation in the central region. Expecta¬ 
tions for the fixed-edge bottom-squark pair-production scenario of Fig. [l] (left), with mass com¬ 
binations o) = (225, 150) GeV, (350, 275) GeV, and (400, 150) GeV for the bottom squark 

and X 2 ' are a ^ so shown. The first scenario presents the illustration of a model that can easily 
be excluded, while the other two present examples of models that are consistent with our data. 



CMS 19.4 fb' 1 (8 TeV) 



m n [GeV] 


Figure 5: Comparison between the observed and estimated SM background dilepton mass 
distributions in the (left) central and (right) forward regions, where the SM backgrounds are 
evaluated from control samples (see text) rather than from a fit. The rightmost bins contain 
the overflow. The vertical dashed lines denote the boundaries of the low-mass, on-Z, and high- 
mass regions. The lower plots show the ratio of the data to the predicted background. The error 
bars for both the main and lower plots include both statistical and systematic uncertainties. 


The results from the dedicated on-Z counting experiment mentioned in Section|4]are presented 
in Tables [ij and [ 5 ] for events with Nj ets > 2 and Nj ets > 3, respectively. The corresponding E™ ISS 
distributions are shown in Fig. [7} The data are seen to agree with the SM predictions across the 
full Ej iss spectrum. 
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Table 3: Results of the edge-search counting experiment for event yields in the signal regions. 
The statistical and systematic uncertainties are added in quadrature, except for the flavor- 
symmetric backgrounds. Low-mass refers to 20 < mu < 70GeV, on-Z to 81 < mu < 101 GeV, 
and high-mass to mu > 120 GeV. 



Low-mass 

On-Z 

High- 

mass 


Central 

Forward 

Central 

Forward 

Central 

Forward 

Observed 

860 

163 

487 

170 

818 

368 

Flavor-symmetric 

722 ± 27 ± 29 

155 ± 13 ± 10 

355 ± 19 ± 14 

131 ±12 ±8 

768 ± 28 ± 31 

430 ± 22 ± 27 

Drell-Yan 

8.2 ± 2.6 

2.5 ±1.0 

116 ±21 

42 ±9 

2.5 ± 0.8 

1.1 ±0.4 

Total estimated 

730 ± 40 

158 ± 16 

471 ± 32 

173 ± 17 

771 ± 42 

431 ± 35 

Observed—estimated 

1302« 

c+20 

°-20 

16 +37 

10 -38 

o+20 

°-21 

47+49 

*'-50 

-62 tg 

Significance 

2.6(7 

0.3 cr 

0.4(7 

<0.1(7 

0.9cr 

<0.1(7 




Figure 6: Data compared with SM simulation for the SF (left) and OF (right) event samples in 
the central region. Example signal scenarios based on the pair production of bottom squarks 
are shown (see text). In the ratio panel below each plot, the error bars on the points show the 
statistical uncertainties of the data and MC samples, while the shaded band indicates the MC 
systematic uncertainty. 
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8 Uncertainties in signal modeling 


CMS 


19.4 fb' 1 (8 TeV) 


CMS 


19.41b' 1 (8 TeV) 




pmiss 


[GeV] 


Figure 7: The E™ 1SS distributions for the on-Z signal regions with (left) >2 jets and (right) 
>3 jets. The uncertainty band shown for the ratio includes both statistical and systematic un¬ 
certainties. Additionally, E™ ISS distributions are drawn for two choices of masses in the GMSB 
scenario. The rightmost bins contain the overflow. 


Table 4: Event yields in the signal region for the dedicated on-Z counting experiment with 
A/jets > 2. Both statistical and systematic uncertainties are included for the background esti¬ 
mates. Signal yields are also shown for two choices of masses (in GeV) in the GMSB scenario 
(statistical uncertainties only). 


E™ 1SS (GeV) 

100-200 

200-300 

>300 

DY background 

336 ± 89 

28.6 ±8.6 

7.7 ±3.6 

FS background 

868 ± 57 

45.9 ± 7.3 

5.1 ±2.3 

Total background 

1204 ± 106 

74.5 ± 11.3 

12.8 ±4.3 

Data 

1187 

65 

7 


GMSB signal yields 


nig = 900, m^o = 150 

22.1 ±0.4 

11.1 ±0.3 

7.2 ± 0.2 

nig = 1100, m^o = 800 

1.1 ±0.04 

1.6 ± 0.05 

7.6 ±0.1 


8 Uncertainties in signal modeling 

Systematic uncertainties associated with the estimation of the SM background are discussed 
in Section [5j This section describes uncertainties in the signal modeling. The impact of the 
uncertainties on the considered signal models is shown in Table [6] 

The uncertainty related to the measurement of the integrated luminosity is 2.6% H36l . The un¬ 
certainty related to the parton distribution functions (PDF) is evaluated using the PDF4LHC 
recommendations 137114111 and amounts to 0-6% in the signal acceptance. Uncertainties in the 
modeling of the reconstruction, identification, and isolation efficiencies amount to 1% per lep¬ 
ton. The uncertainty of the corrections to account for lepton reconstruction differences between 
the fast and full simulations is 1% per lepton. The uncertainty in the trigger efficiencies is about 
5%. Uncertainties in the muon momentum scale are negligible, whereas for electrons the un¬ 
certainty in the energy scale is 0.6% (1.5%) in the barrel (endcap) region, leading to a 0-5% 
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Table 5: Event yields in the signal region for the dedicated on-Z counting experiment with 
A/jets > 3. Both statistical and systematic uncertainties are included for the background esti¬ 
mates. Signal yields are also shown for two choices of masses (in GeV) in the GMSB (statistical 
uncertainties only). 


E™ 1SS (GeV) 

100-200 

200-300 

>300 

DY background 

124 ± 33 

12.7 ±3.8 

3.2 ±1.8 

FS background 

354 ± 28 

26.5 ± 5.4 

2.0 ±1.4 

Total background 

478 ± 43 

39.2 ± 6.6 

5.3 ± 2.3 

Data 

490 

35 

6 


GMSB signal yields 


trig = 900, m^o = 150 

22.0 ± 0.4 

11.0 ±0.3 

7.1 ±0.2 

rtig = 1100, m^o = 800 

1.1 ± 0.04 

1.5 ±0.05 

7.4 ± 0.1 


Table 6: Summary of systematic uncertainties for the signal efficiency. 


Uncertainty source 

Impact on signal yield [%] 

Luminosity 

2.6 

PDFs on acceptance 

0-6 

Lepton identification/isolation 

2 

Fast simulation lepton identification/isolation 

2 

Dilepton trigger 

5 

Lepton energy scale 

0-5 

£miss 

0-8 

Jet energy scale/resolution 

0-8 

ISR modeling 

0-14 

Additional interactions 

1 


uncertainty in the signal yield. The uncertainty in the jet energy scale l42l varies between 0- 
8 %. The uncertainty in the jet energy scale is propagated to evaluate the uncertainty associated 
with the Ej' ss distribution, which is found to be 0-8%. The uncertainty associated with the 
modeling of initial-state radiation (ISR) Il43l is 0-14%. The uncertainty in the correction to ac¬ 
count for pileup in the simulation is evaluated by shifting the inelastic cross section by ±5%. 
The impact on the signal yield is found to be about 1%. 


9 Interpretation 

Based on the results of the counting experiments, exclusion limits are determined. The kine¬ 
matic fit is not used for this purpose because it assumes a specific shape for the dilepton mass 
spectrum of signal events and so is more model dependent. The limits are calculated at 95% 
confidence level (CL) using the CL S criterion with the LHC-style test statistic Il34l j44', l45l . tak¬ 
ing into account the statistical and systematic uncertainties in the signal yields and the back¬ 
ground predictions discussed in Sections [8] and |5} respectively. The different systematic un¬ 
certainties are considered to be uncorrelated with each other, but fully correlated among the 
different signal regions. The systematic uncertainties are treated as nuisance parameters and 
are parametrized with log-normal distributions. 

For the fixed- and slepton-edge scenarios all six signal regions (low-mass, on-Z, and high- 
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10 Summary 


mass, for the central and forward lepton regions) are combined. The resulting exclusion limits 
for these scenarios are shown in Fig. [8] The production of bottom squarks is considered, as 
the excess observed in data consists predominantly of events with at least one b jet. In the 
fixed-edge scenario, the mass difference between the %% an d T? neutralinos is fixed to 70 GeV, 
resulting in an edge in the nifj, spectrum at this value. This choice for the edge position is 
motivated by the observed excess in the low-mass region of the counting experiment. Bottom- 
squark masses between 200 and 350 GeV are probed, depending on the value of the mass. 
In the slepton-edge scenario, the mass is set to 100 GeV, leaving the position of the edge as a 
free parameter that approximately corresponds to the mass difference between the and Xy 
The branching fraction into dilepton final states is larger than in the fixed-edge scenario, and 
bottom-squark masses between 450 and 600 GeV are probed, again depending on the value of 
the X 2 mass. The loss of sensitivity seen in Fig. [8] (right) for X 2 masses around 250 GeV occurs 
because the peak of the triangular signal shape lies within, or close to, the gaps in acceptance 
between the low-mass and on-Z, or the on-Z and high-mass signal regions. 
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> 

CD 


0 


o~c\j 
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1 
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_Q 

Cl 
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CMS 


700 F 
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600- 


400 
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Figure 8: Exclusion limits at 95% CL for the fixed- (left) and slepton-edge (right) scenarios in 
the mg-m ~0 plane. The color indicates the excluded cross section for each considered point 
in parameter space. The intersections of the theoretical cross section with the expected and 
observed limits are indicated by the solid and hatched lines. The 1 standard deviation (a) 
experimental and theoretical uncertainty contours are shown as dotted lines. 


The results from the dedicated on-Z signal regions are used to place limits on the GMSB sce¬ 
nario. In this scenario, there are two free parameters: the masses of the gluino (mg) and the 
As signal events typically have large jet multiplicities, the exclusive bins requir¬ 
ing Njets > 3 and E™ 1SS in the ranges 100 < E™ ss < 200 GeV, 200 < Ey ISS < 300 GeV, and 
Ey ISS > 300 GeV are used. The results are shown in Fig. 9jin the plane of ni-o versus mg. These 
results probe gluino masses up to about 900-1100 GeV depending on the x^ mass. The limit is 
least stringent when m ~0 is close to the Z boson mass. 


10 Summary 

We have presented a search for physics beyond the standard model in the opposite-sign dilep¬ 
ton final state using a data sample of pp collisions collected at a center-of-mass energy of 8 TeV, 
corresponding to an integrated luminosity of 19.4 fb -1 , recorded with the CMS detector in 
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Figure 9: Exclusion limits at 95% CL for the GMSB scenario in the ntg-m^o plane. The color 
indicates the excluded cross section for each considered point in parameter space. The inter¬ 
sections of the theoretical cross section with the expected and observed limits are indicated by 
the solid and hatched lines. The 1 standard deviation (n) experimental and theoretical uncer¬ 
tainty contours are shown as dotted lines. 


2012. Searches are performed for signals that either produce a kinematic edge, or a peak at the 
Z boson mass, in the dilepton invariant mass distribution. For regions dominated by flavor- 
symmetric backgrounds, i.e., backgrounds that produce opposite-flavor events (e' y e }i ' ) 
as often as same-flavor events (e - e , ji~}i "), we achieve a precision of about 5% (10%) for the 
estimated number of standard model background events in the central (forward) lepton rapid¬ 
ity regions. We do not observe evidence for a statistically significant signal. The maximum 
deviation from the null hypothesis is at the level of 2.6 standard deviations and is observed in 
the dilepton mass window 20 < mu < 70 GeV. 

We interpret the results of the search for a kinematic edge in the context of simplified models 
consisting of bottom-squark pair production, with each bottom-squark b decaying to a bottom 
quark and the y, neutralino. Exclusion limits are set in the in^-m-o mass plane for two scenar¬ 
ios. In the fixed-edge scenario, the mass difference between the and y)’ particles is fixed to 
70 GeV. Bottom-squark masses between 200 and 350 GeV are excluded at 95% confidence level 
in this scenario. In the slepton-edge scenario, the y!) decays to an on- or off-shell Z boson and 
the y ( | lightest supersymmetric particle or to a slepton and a lepton, with a 50% probability for 
each possibility. Bottom-squark masses between 450 and 600 GeV are excluded at 95% confi¬ 
dence level in this second scenario. In both scenarios, the sensitivity depends on the mass of 
the X 2 ■ Finally, a dedicated search for events containing an on-shell Z boson is interpreted in a 
model of gauge-mediated supersymmetry breaking, in which the Z bosons are produced in de¬ 
cay chains initiated through gluino pair production. Gluino masses between 900 and 1100 GeV 
are excluded at 95% confidence level, depending on the mass of the lightest neutralino y. 
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